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Abstract 

Based on flat Friedmann- Robert son Walker cold dark matter (CDM) type models driven 
by non-relativistic matter and an exotic fluid (quintessence) with an equation of state: pq = 
wpQ (—1 < w < 0), we investigate whether or not the large scale dynamical effects regarding 
the cluster formation and virialization are related to the cosmic equation of state. Using the 
non-linear spherical collapse we find that the cluster formation rate, in quintessence models, 
is intermediate between the open and ACDM respectively. For the QCDM case, using the 
virial theorem and energy conservation and assuming a spherical mass overdensity shell, we 
obtain analytically the ratio between the final (virial) and the turn-around radius. We find 
that the above ratio is almost independent from the equation of state. 
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1. Introduction 

Recent advances in observational cosmology have strongly indicated that we are living 
in a flat, accelerating universe with low matter (baryonic and dark matter) density (Riess 
et al. 1998; Perlmutter et al. 1999; de Bernadis et al. 2000; Efstathiou et al. 2002; Persival 
et al. 2002 and references therein). The available high quality cosmological data (Type la 
supernovae, CMB, etc.) are well fitted by an emerging "standard model" , which contains cold 
dark matter (CDM) to explain clustering and an extra component with negative pressure 
(dark energy), usually named "quintessence", to explain the inflationary flatness prediction 
(f^tot = 1) as well. 

The last few years there have been many theoretical speculations regarding the nature 
of the above exotic dark energy. Most of the authors claim that a time varying A-parameter 
(cf. Ozer & Taha 1987), a scalar field which rolls down the potential V(<j)) (Ratra & Peebles 
1988; Frieman et al. 1995 and references therein) or an extra "matter" component, which is 
described by an equation of state pq = wpQ where — 1 < w < 0, could be possible candidates 
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for quintessence. As a particular case, A-models can be obtained from quintessence models 
for w — — 1. An excellent paper by Tsagas (2001) gave a different perspective to the problem 
claiming that maybe the effects of magnetic fields could resemble that of dark energy through 
the coupling between magnetism and curvature of spacetime. While a variety of observations 
indicate that w < —0.6 for a flat universe (cf. Efstathiou 1999; Dev, Sethi, & Lohiya 2001; 
Sereno 2002; Ettori, Tozzi & Rosati 2002 and references therein). In this work, we investigate 
the cluster formation processes utilizing the non-linear spherical model for a family of various 
cosmological models in order to understand better the theoretical expectations of negative 
pressure models as well as the variants from the ACDM case. 



2. Basic Equations 

For homogeneous and isotropic cosmologies, driven by non relativistic matter and an 
exotic fluid (quintessence models) with equation of state, pq = wpQ with — 1 < w < 0, the 
Einstein field equations can be given by: 

a\ 2 8tiG k 

(Pm + PQ)--z (!) 



a 
and 

= -AnG[(w + \)p Q + \p m \ , (2) 

where a(t) is the scale factor, p m oc (1 + z) 3 is the matter density and pq oc (1 + z) 13 is the 
dark energy density [f3 = 3(1 + w)\. While k = —1, or 1 for open, flat and closed universe 
respectively. Thus, the scale factor evolves according to Friedmann equation: H 2 = (a /a) 2 . 
In order to transform the latter equation from time to redshift we utilize the following 
expression: 

— = 1 ( 3 ) 

dz H E(z)(l + z) ' V ; 

with 

e(z) = [n m (i + zf + n k (i + zf + n Qo (i + z f] 1/2 (4) 

where the Hubble parameter is given by: H(z) = H Q E(z), while Q m = 8irGp /3H 2 (den- 
sity parameter), Q k = —k/H a (curvature parameter), Qq = 8nGpQ /3H 2 (dark energy 
parameter) at the present time, which satisfy f2 m + Qk + Qqo = 1 and H Q is the Hubble 
constant. 

In addition, to £l m (z) also Qq(z) could evolve with redshift as 

^(^) = E2 ^ and n Q (z) = . (5) 
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Note, that A-models can be described by quintessence models with w strictly equal to -1 
while if w = the equation of state behaves like that of pressureless matter. It is interesting 
to mention that different values of w could yield flat = 0) cosmological models for which 
there is not a one-to-one correspondence between the global geometry and the expansion of 
the universe. Indeed, in flat low-f2 m with w = —1/3 model, the functional form of the dark 
energy density can be given by pq oc (1 + z) 2 . Therefore, the equation of state pq = wpQ 
leads to the same expansion as in an open universe. In other words, pq plays a similar role 
to the curvature, despite the fact that this quintessence model has a spatially flat geometry! 
Also in the case of a flat low-fi m with —1 < w < —1/3 orO < (3 < 2 model, we have 
Pq oc a" 13 < a~ 3 oc p m , which means that the dark energy density falls off at a slower rate 
than cold dark matter. This is very important because the dark energy component starts 
to dominate the mass density in the universe, especially at the late times, thus creating an 
accelerating expansion. 

As for the power spectrum of our CDM models, we consider P(k) ~ k n T 2 (k) with 
scale-invariant (n = 1) primeval inflationary fluctuations. In particular, we consider three 
different spatially flat low-f2 m cosmological models with negative pressure, where: (1) an 
Vl m = 1 - n A = 0.3 with w = -1 (ACDM), (2) a model with fi m = 1 - n Qo = 0.3 
with w = -2/3 (QCDM1) and to this end (3) an fi m = 1 - Q Qo = 0.3 with w = -1/3 
(QCDM2). Our models have T ~ 0.2, in approximate agreement with the shape parameter 
estimated from galaxy surveys (cf. Maddox et al. 1990) and Hubble constant H a = 100/i 
km s _1 Mpc _1 (where h = 0.65). Finally, the latter cosmological models are normalized to 
have fluctuation amplitude in 8 /i -1 Mpc scale of a 8 = 0.50 ± 0.1f2~ 7 (Wang & Steinhardt 
1998) with 7 = 0.21 - 0.22w + 0.33fi m . 



The growth factor as a function of redshift 1 for the mass density contrast, 5 = (5p m /p m ), 
modeled as a pressureless fluid (cf. Peebles 1993; Lokas 2001 and references therein) is: 



In case of w — —2/3, D(z) can be obtained by Wang & Steinhardt (1998; their equation 
14). Below we attempt to investigate how the growing mode of perturbations is affected by 
the time evolution of the dark energy contribution to O to t ■ To do so we follow the ideas of 



1 D(z) = (1 + z) 1 for an Einstein-de Sitter universe. 



3. Cluster formation in QCDM Cosmologies 




(6) 
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Lightman & Schechter 1990, Lahav et al. (1991), and Carrol, Press & Turner (1992), who 
use a convenient approximation for D(z): 



D(z)= 5fim 



Q^-Q Q + (l + ^Q m )(l + AQ Q ) 



(7) 



2(1 + z) 

where we leave unconstraint the values of a and A and we use the generic expression for a, 
defined by the Wang & Steinhardt (1998): 

5-w/(l-w) 125 (l-6u>/5) 3 V ; KJ 

Utilizing a x 2 minimization procedure between the complete D(z) solution and equations 
(7), (8), considering different values for A and w for fl m = 0.3, we obtain a roughly quadratic 
relation for A, fitted by: A ~ 1.742 + 3.343w + 1.615w 2 . Note that for the ACDM case 
(w = —1), the result A ~ 0.014 with a ~ 6/11 agrees with Lahav et al. (1991) who 
found A ~ 1/70 and a ~ 0.6. Finally, for the quintessence models QCDM1 (w = —2/3) 
and QCDM2 (w = -1/3) we have (A, a) ~ (0.232,0.565) and (A, a) ~ (0.810,0.584) 
respectively. In Figure 1, lines represent the complete growing mode solution in contrast 
to the D(z) approximation formula, which is represented by the points. Thus, it becomes 
evident that the growing mode approximation works extremely well. 

The concept of estimating the fractional rate of cluster formation, considering the non- 
linear spherical collapse model, has been brought up by different authors (cf. Peebles 1984; 
Weinberg 1987; Martel & Wasserman 1990; Lahav et al. 1991; Richstone, Loeb & Turner 
1992). Here, we present the basic steps of the cluster formation processes that we will use, fol- 
lowing the above ideas. Basically, the above authors introduced a procedure which computes 
the rate at which mass join virialized structures, which grow from small initial perturbations 
in the universe. It is well known that the basic cosmological equations, mentioned before, 
are correct either for the entire universe or for homogeneous spherical perturbations [by re- 
placing the scale factor with radius R(t)]. Therefore, having bound perturbations which do 
not expand forever the time they need to recollapse (at some redshift z = Zf) is twice the 
turn-around time if = 2t ta . 

Furthermore, making the assumption that when the matter epoch just begins, the uni- 
verse is described by (i) an unperturbed Hubble flow, (ii) a matter fluctuation field which 
has a Gaussian distribution: 

^ (A) = v^ exp (^) <iA ' (9) 

where A is the mass density contrast and it contains a dark energy component. The rms 
mass fluctuation amplitude at 8 /i _1 Mpc can be expressed as a function of redshift as a(z) = 
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D(z)crg. Then, for the bound perturbations that we consider, following the standard cluster 
formation pattern (cf. Richstone et al. 1992), we can obtain the ratio of the collapse time 
to the current age of the (unperturbed) universe: 



t _ f 00 dx 
T~ ~ 2 J Z (l+x)M(A,x) 



dx 



-i -i 



o (l + x)E(x) 



(10) 



where 



M(A, z) = [Q m (l + A)(l + zf + (1 - n m (l + A) — 

n Qo )(i + z ) 2 + n Qo (i + zf} 1 / 2 . 

The solution of the integral of eq.(9), which is Gaussian, describes the fraction of the 
universe (characterized by f2 m , Qq q and cr 8 ) on some specific mass scale that has already 
collapsed at time t and is given by (see also Richstone et al. 1992): 



1 - erf 



M(t/T ) V 
V V2a )_ 



(11) 



The next step is to normalize the probability to give the number of clusters which have 
already collapsed by the epoch t (cumulative distribution), divided by the number of clusters 
which have collapsed at the present epoch: 



r = wit h 



/ = 



(n)M 
p c f2 m 



(12) 



Where p c ~ 2.78 x 10 u h 2 M Q Mpc -3 is the critical density, (n) is the number density of 
clusters which have collapsed prior to the present epoch. The parameter (n) can be defined 
utilizing the Abell/ACO cluster catalog, which is a volume-limited sample within ~ 180 — 
200/i" 1 Mpc and which a general agreed value is: (n) ~ 1.8 x I0~ 5 h 3 Mpc~ 3 . Therefore, 
considering virialized clusters of the mass scale of rich Abell clusters, M ~ 1O 15 /i -1 M , it is 
a routine to obtain the ratio of collapsed matter at the present time /(1O 15 M ) = O.OGSfi" 1 . 

It is obvious that the above generic of form eq.(12) depends on the choice of the back- 
ground cosmology. Indeed the relationship between A, z, t/T Q and as is different in different 
cosmologies (cf. Mo & White 1996; Magliocchetti et al. 2001), where: 



A(z) = 6 C _5c_ 
a 8 D(z)a 8 a{z) 



(13) 



The value 5 C = 1.686 corresponds to the spherical top-hat model in f2 m = 1, but it has been 
shown that S c depends only weakly on the cosmology (Eke, Cole & Frenk 1996). Considering 
the three low-f2 m spatially flat models, described in section 2 we can obtain complementary 
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predictions also for the old "standard" CDM model with Q m — 1, h — 0.5 (SCDM) and the 
open model with Q m = 0.3 and h = 0.65 (OCDM). The latter two cosmological models are 
normalized by the observed cluster abundance at zero redshift; a 8 = 0.55f2~ 0,6 (Eke et al. 
1996). 

In figure 2 we present the behavior of eq.(12) as a function of redshift and t/T . The 
behavior of the cluster formation rate, J 7 , has the expected form, ie. it is a decreasing 
function of redshift. For the high density universe Q m = 1 we found the known behavior 
(Richstone et al. 1992) in which galaxy clusters started to form only very recently while in 
an open or a flat low-density universe, clusters should appear to be formed at z ~ 2 due to 
the fact that clustering effectively freezes at high redshifts. In this framework, we have found 
that the cluster formation rate in quintessence models has an intermediate growth between 
that of an open and ACDM models respectively. Finally, it is obvious that for w — ► —1 the 
quintessence cluster formation pattern tends to that of ACDM case, as it should. 



3.1. Cluster virialization with dark energy 

In this section, based on the notations of Lahav et al. (1991), Wand & Steinhardt (1998) 
and Lokas & Hoffman (2003) we study the cluster virialization in QCDM models. Here we 
review only some basic concepts of the problem. Assuming a spherical mass overdensity 
shell, utilizing both the virial theorem T = —\Ug + Uq and the energy conservation Tf + 
Ugj + Uqj = Ua,ta + Uqjo, where, T is the kinetic energy, Uq is the potential energy and 
Uq is the potential energy associated with quintessence. Therefore, we can obtain a cubic 
equation which relates the ratio between the final (virial) i?f and the turn-around outer 
radius R t &. 



2nu \Rt) ~ (2 + nt) V^W +1 = ' (14) 

with an approximate solution (see Wand & Steinhardt 1998): 

Rt_ 1 - n u /2 
i? ta 2 + n t - 3n u /2 

where n u = 2fi Qo (l+2; f ) 3 ( 1+ "')/Cfi m (l+^a) 3 , n t = 2{l Qo (l+z tll ) 3w /({l m and C = p d (*ta)A»m(*ta 
and for — 1 < w < is a weakly f2 m -dependent (see Wand & Steinhardt 1998). In this study, 
we derive analytically the exact solution of the above cubic equation, having polynomial 
parameters: ai = 0, a 2 = —(2 + n t )/2n u and a 3 = l/2n u . Indeed, using some basic elements 
from the Algebra (see the Appendix) the discriminant of eq.(14) is: 

2(2 + n t f - 27n u 



V(n u ,n t ) = 
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Owing to the fact that we are living in an accelerating universe £Iq > with low matter 
(baryonic and dark matter) density, the condition for an overdensity shell to turn around is 
< n t < n u < 1 which gives V{n u ,n t ) > and therefore all roots of the cubic equation are 
real (irreducible case) but one of them 

r 2 = -^- = ~[i2 + V3M] (15) 

corresponds to expanding shells (for R, M and the other parameters see the Appendix). 
Thus, according to the solution described above, figure 3 shows the surface behavior of the 
exact virial solution in parametric form with: 

i; , / 27 2 27£>\ 1/6 , 6 1 (2n u y/ZD\ 
r = {W u + —) and 3 = 3 arCtan {—^- ) 

It is interesting to say that the minimal possible ratio is Rf/R ta — 0.35 which is similar to 
the ACDM solution derived by Lahav et al. (1991) for which case we have n u = n t = n = 
f2A/C^m(l + <£ta) 3 ) where £7a is the cosmological constant parameter at the present time. 
The latter minimal ratio implies that the final stage of the cluster virialization is almost 
independent from the equation of state. 



4. Conclusions 

Using the non-linear spherical collapse model, we derive the cluster formation rate as a 
function of redshift in CDM models with negative pressure. We verify that the formation rate 
in quintessence models ( — 1 < w < 0) is an intermediate case between the open and ACDM 
models respectively. Furthermore, we study the cluster virialization in QCDM models and 
we derive the exact solution. The minimal possible ratio is R{/R ta , ~ 0.35 which is consistent 
with the ACDM solution derived by Lahav et al. (1991) which means that the above ratio 
is almost independent from the equation of state. 



Appendix 

Without wanting to appear too pedagogical, we remind the reader of some basic elements 
of Algebra. Given a cubic equtation: x 3 + a^x 1 + a 2 x + a 3 = 0. Let D be the discriminant: 

V = a\a\ - Aa\ - 4afa 3 - 27a 2 + 18aia 2 a 3 (16) 

and 

3 9 27 3V3V 
xi = -a ± + -a x a 2 - —a z , x 2 = — 
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then, the roots of the above equation are: 



ri = -^ + i[?i + ft] (17) 



r2 = -j + l[e 2 qi + eq 2 ] (18) 

r 3 = -j + \[eqi + e 2 q 2 ] (19) 

where q n = {x\ ± ixz) 1 ^ and e = ~ 1+ 2 ^~~^ ■ If X> < 0, we have one real root {r\) and a pair of 
complex conjugate roots. If T> — 0, all roots are real and at least two of them are equal. If 
V > 0, all roots are real (irreducible case). In that case r\, r 2 and r 3 can be written: 

n = -| + f* (20) 



r 2 = -^-^[i2 + V3M] (21) 



r 3 = -|-i[ J R-v / 3M] (22) 



where i? and M are given by 



i? = r 1/3 cos(-) and M = r 1/3 sin(-) (23) 
3 3 



with r = \J x\ + x\ and 6 = Arctan(x 2 /a;i). 
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Fig. 1. — Growth factor for matter density fluctuations as a function of redshift in linear the- 
ory using cosmological models (Q m = 0.3) with negative pressure. Note that lines represent 
the complete D(z) solution while points represent the D(z) approximation formula. 
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Fig. 2. — Theoretical predictions of the fractional rate of cluster formation as a function of 
redshift and fractional time t/T a , for different cosmo logical models. 
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Fig. 3. — The ratio of the final (virial) to turn-around radius of a virialized cluster utiliz- 
ing a QCDM cosmology. We give the surface solution in parametric form (r 1//3 , 6/3 see the 
Appendix). 



